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Abstract. Following the use of approximate symmetries for the Schwarzschild spacetime 
by A.H. Kara, F.M. Mahomed and A. Qadir (Nonlinear Dynam., to appear), we have in- 
vestigated the exact and approximate symmetries of the system of geodesic equations for 
the Reissner-Nordstrom spacetime (RN). For this purpose we are forced to use second order 



approximate symmetries. It is shown that in the second-order approximation, energy must 
be rescaled for the RN metric. The implications of this rescaling are discussed. 
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1 Introduction 

r> . 

Though most of the well known spacetimes of General Relativity (GR) are time-translational 
invariant, generally this is not guaranteed [2]. Consequently energy is not generally conserved 
and hence there is no good definition for energy [3]. Nowhere is this problem faced as severely 
as in gravitational wave spacetimes, where energy is not given by the stress-energy tensor, which 
is zero [2] . For these spacetimes one can check that they impart energy to test particles in their 
path [U [5] despite the fact that the stress-energy tensor is zero. The lack of a good definition of 
energy, also leads to problems with the definition of mass [3]. As such, one needs to use some 
concept of time symmetry that allows for slight deviations away from exact symmetry. 

This approach was attempted earlier by various people. There have been a number of different 
definitions of "approximate symmetry". One idea was to assume that conservation of energy 
holds asymptotically [6] and to examine whether it would work for gravitational radiation and to 
define a positive definite energy. This seems unsatisfactory as the gravitational energy should 
reach (the place) infinity. There may then be problems with orders of approximation being 
consistent. An altogether different approach was taken by providing a measure of the ex- 
tent of break-down of symmetry. The integral of the square of the symmetrized derivative of 
a vector field was divided by its mean square norm [8]. This led to what was called an 
almost symmetric space and the corresponding vector field an almost Killing vector [9]. This 



*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
http://www.emis.de/journals/SIGMA/symmetry2007.html 



2 



I. Hussain, F.M. Mahomed and A. Qadir 



measure of "non-symmetry" in a given direction was applied to the Taub cosmological solu- 
tion [TU] and to study gravitational radiation. It provides a choice of gauge that makes calcu- 
lations simpler and was used for this purpose [llj . Essentially based on the almost symmetry, 
the concept of an "approximate symmetry group" was presented [12]. However it has not been 
unequivocally successful either. The approach of a slightly broken symmetry seems promising, 
but merely providing simplicity of calculations is not physically convincing. Other approaches 
need to be tried, to find one that seems significantly better than others. 

The development of geometry has been driven by its application to kinematics and dyna- 
mics. As such, one might look at a geometrically driven definition of symmetry. On the other 
hand, we are concerned with conservation laws, which are given by the invariants of the Euler- 
Lagrange equations. Hence, for our purposes, Lie symmetries as embodied in isometries and 
in Noether's theorem, should provide us the desired approach to define "approximate sym- 
metries". Now, there is a connection between isometries and the symmetries of differential 
equations (DEs) through the geodesic equations [HKllj. We propose that the Baikov-Gazizov- 
Ibragimov concept of "approximate symmetry" [151 [T6] could be extended and adapted for the 
purpose of defining energy in GR via "approximate isometries" through the above-mentioned 
connection between the symmetries of geometry and the symmetries of differential equations. 
Using that connection approximate symmetries of the Schwarzschild metric were discussed [1] . 
The Schwarzschild metric has much fewer symmetries than the Minkowski metric (only 4 gen- 
erators). One would expect that in the limit of small gravitational mass the "lost" symmetries 
should be "approximately" recovered. Of course, this is a static metric and energy is conserved 
and it is not relevant for a discussion of energy. However, linear and spin angular momentum 
conservation is lost. The "trivial" approximate symmetries of the geodesic equations recover 
these as approximately conserved quantities. To be able to interpret the approximate symmetry 
results so obtained we need to first follow how higher order approximate symmetries are to be 
understood for our purpose and then how the lack of time-translational symmetry on a globally 
defined spacelike hypersurface is to be interpreted. As such, before dealing with gravitational 
waves it is necessary to study the Reissner-Nordstrom (RN) metric and the Kerr metric from 
this point of view. In this paper we present the discussion of RN spacetime leaving the other 
spacetimes for subsequent consideration. 

In going to the RN spacetime the approximate symmetries and conservation laws, which were 
recovered in the first-order approximation [1], are lost. Nor is there any approximate symmetry 
for the reduced orbital equation. However one would expect that in the limit of small charge we 
should recover the lost symmetries. To recover the conservation laws lost again we need to appeal 
to second-order approximate symmetries. It makes no difference whether the exact or perturbed 
equation is used in the definition of first-order approximate symmetries. In this paper it is shown 
that in the definition of second-order approximate symmetries it makes a significant difference. 
We use the perturbed equation in the definition of second-order approximate symmetries, the 
expression for energy of a test particle in the RN metric is re-scaled. 

The plan of the paper is as follows. In the next section we explain the concept of symmetry 
and approximate symmetry of DEs, and approximate symmetries of the Schwarzschild metric. 
In Section 3 approximate symmetries of the orbital equation for the RN metric and of the system 
of geodesic equations for this metric are discussed. Finally a summary and discussion is given 
in Section 4. 

2 Symmetries and approximate symmetries of DEs 

Symmetries are very useful because they are directly connected to the conservation laws through 
Noether's theorem {17} IT8| I19j . If, for a given system of DEs there is a variational principle, 
then a continuous symmetry under which the action functional remains invariant yields a con- 
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servation law [20\ [2~T| [22] [23] . Here we are interested in Lie symmetry methods [24^ . The set 
of infinitesimal symmetry generators of a system of differential equations form a Lie algebra 
[2 H 1251 [16] . Geometrically the symmetries of a manifold are characterized by its Killing vectors, 
or isometries, which also form a Lie algebra [26]. The symmetries of the manifold are inherited 
by the geodesic equations on it with additional symmetries [13] • They give quantities conserved 
under geodesic motion [27] and lead to first integrals of the geodesic equations [28] . 

It may often happen that a manifold does not possess exact symmetry but approximately does 
so. We may be able to obtain more information from near symmetry or broken symmetry, than 
from the exact symmetry always maintained. As such approximate symmetries of a manifold 
are worth exploring. Methods for obtaining the approximate symmetries of a system of ordinary 
differential equations (ODEs) are available in the literature [29} 115]. 

A system of p ODEs of order n \25\ 116] 

E( S ;x( S ),x'( S ) 1 x"( S ),..,x(»)( S ))=0, 
has the symmetry generator 

X [n] =^)§- S + V a (s,x)^+V? s (s,^ (1) 

if X[ n l annihilates the system of equations E (for solutions of the equation) 
XN (E)| E=0 = 0, 

where x is a point in the underlying n-dimensional space and x(") is the nth derivative with 
respect to s. The prolongation coefficients are given by 

a _ a(n) > -i wVipt-p n a — n a 

V,(n)- ds x dg , n>L, wnere rj^-r]^ 

If p = 1 then the system reduces to a single equation. 

The fcth-order approximate symmetry of a perturbed system of ODEs 

E = E + eEi + e 2 E 2 + • • • + e k B k + 0(e k+1 ) (2) 
is given by the generator 

X = X + eXi +e 2 X 2 + • • • + e k X k , (3) 
if the symmetry condition 

XE := [(X + eXi + e 2 X 2 + • • • + e fc X fc )(E + eE 1 + e 2 E 2 + • • • + e fc E fe )]|2} = 0(e k+1 ) (4) 

is satisfied [191 [23 EQ] . Here Eo is the exact system of equations, Ei is the first-order approx- 
imate part and E 2 is the second-order approximate part of the perturbed system and so on; 
Xo is the exact symmetry generator, Xi the first-order approximate part, X 2 the second-order 
approximate part of the symmetry generator and so on. These approximate symmetries do not 
necessarily form a Lie algebra but do form a so-called "approximate Lie algebra" |31j . 

If ([2} admits ([3]) with Xo 7^ 0, then Xo is an exact symmetry of the unperturbed Eo = 
and Xo is a stable symmetry for a given perturbation. The perturbed equation ([2]) always has 
the approximate symmetry eXo which is known as a trivial symmetry and X given by ([3]) with 
Xo 7^ and Xi 7^ Xo is called a non-trivial approximate symmetry [32J. 

An alternate method for defining approximate symmetries of DEs was given by Fushchich 
and Shtelen [33J. They interchange the order of approximation and take the limit, between 
the parameter of the symmetry generator of the algebra on one hand and the approximation 
parameter on the other hand. This method is compared with that of Baikov et al. in |34[ l35j. For 
our purpose it makes no difference which of these methods is used. We will follow the method 
of Baikov et al. 1151. 
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2.1 Symmetries and approximate symmetries of the Schwarzschild metric 

Minkowski spacetime has the Poincare symmetry algebra so(l,3) © S K 4 (where S denotes semi 
direct sum) with the 10 generators \36\ [2] 

X = ^, Xi = cos^-cat0sin^, (5) 

X 2 = sin ^+ cot0 cos <j>^-, x 3 = ^, (6) 

d cos cos 4> 6 esc sin (p 6 

X4 = sm 9 cos cp— H — — , (7) 

or r 69 r oq> 

d cos 9 sin 6 6 esc 9 cos <fi d 

X 5 = sin sin 0— H — H — , (8) 

ar r a0 



r 



sin0 6 

r sin 9 cos 9 / . ,0 cos cos ^ 9 esc sin (f) d \ 

X 7 = — + ci sm0cos<£— + ^ ^-7 , (10) 

c ot \ ar r 09 r 6<p J 

r sin sin <i 9 / . „ ,0 cos sin <A 9 esc cos 4> d \ , , 

X 8 = ^ — + rf sm0sm«/» — + ^— + ^— , 11 

c OT \ ar r 69 r 6tp J 

rcos0 a / „ a sin0 6 \ 

^ = —at +c, \ cosl> ¥r-—m)- (12) 

Xo, X4, X5 and ~Kq are the spacetime translations (with a Lie sub-algebra (M 4 )) which provide 
the laws of conservation of energy and linear momentum, Xi, X2 and X3 are the rotations 
(with a Lie sub-algebra (so(3))) which provide the laws of conservation of angular momentum 
and X7, X§ and X9 are the Lorentz transformations (with a Lie sub-algebra (so(3))) which 
provide the laws of conservation of spin angular momentum; via Noether's theorem |17| . The 
symmetry algebra for the geodesic equations is sZ(6,R), which has many symmetries that do 
not correspond to conservation laws |14| . 

For the Schwarzschild metric, corresponding to a point mass, M, situated at the origin only 
the first 4 generators (Xo, Xi, X2, X3) apply, yielding only conservation of energy and angular 
momentum. The generators X5,...,Xg, yielding conservation of linear momentum and spin 
angular momentum are "lost". Using the approximate symmetries for the system of geodesic 
equations for this metric with e = 2GM/c 2 [l] (where G is Newton's gravitational constant 
and c is the speed of light in vacuum) the "lost" symmetries are recovered as trivial approximate 
symmetries yielding approximate conservation of those quantities. The symmetry algebra of the 
geodesic equations for the Schwarzschild metric is so(3)©M®d2 (where d,2 is the dilation algebra 
with generators 6/ds and sd/6s) [H]. The geodesic equations can be reduced to a single orbital 
equation that has the two non-trivial approximate symmetries given by the generators 

6(6 6 
X a i = sm cp— + e 2 sm <j)— + u cos 4>— 
6u V °4> ou, 

d / ,8 



X a2 = cos <j>- e 2 cos 4>— - u sin T 

6u \ 6<p 6u 

3 Symmetries and approximate symmetries 
of the Reissner— Nordstrom metric 

The field of a point massive electric charge at rest at the origin is given by the RN metric (see 
for example [37] ) 

ds 2 = e u M 2 _ e ~v dr 2 _ r 2 dQ 2 _ r 2 ^2 Q ^2 



Second-Order Approximate Symmetries of the Geodesic Equations 



5 



with 

e* = l-^ + ^, (14) 
c 2 r c 4 r 2 

where Q is the electric charge of the point gravitational source. Electromagnetism is the only 
long range force in Nature other than gravity and this is the only spherically symmetric, static 
exact solution of the "sourceless" Einstein-Maxwell equations. In the chargeless case (Q = 0) 
it reduces to the Schwarzschild metric. It is of interest to look at the symmetry structure of 
this metric and the corresponding symmetries, and approximate symmetries of the geodesic 
equations. 

It had been pointed out pQ that there is a difference between the conservation laws obtained 
for the system of geodesic equations and the single, orbital equation for the Schwarzschild 
metric. It was further remarked that it should be checked if this difference also holds for other 
spacetimes. We investigate this question for the orbital equation in the RN metric, 

f^ + u - 9J}£ _ 9@L U + 3GM U 2 _ 2G Q 2 U 3 

d(p 2 h 2 c 2 h 2 c 2 c 2 

where h is the classical angular momentum per unit mass and u = 1/r. In the classical limit 
c — > oo it gives the classical orbital equation and for Q = it yields the Schwarzschild orbital 
equation. For the approximate symmetries of this equation we take 

2GM GQ 2 2 
6 = -^-' — = ke - (15) 

For an RN black hole (see for example [38]) that is 

M 2 > Q 2 , we have < k < 1/4. 

If k > 1/4 the metric represents a naked singularity. Taking e = in the orbital equation, the 
exact symmetry generators are 



i—, Xi = cos0— , X 2 = sm0— , X 3 = — 
ou ou ou am 



d d d d 

X4 = cos2<p— — u sin 20— , X5 = sin 20— + u cos 20— , 

ocp ou d(p ou 

d d d d 

X 6 = it cos 0— — u 2 sin0— , X7 = usin0— + u 2 cos 0—. 

0(p ou ocp Ou 

Retaining terms of first-order in e and neglecting 0(e 2 ), the first approximate symmetry genera- 
tors are given by (|13p . In the second approximation, that is when we retain terms quadratic 
in e, this equation possesses no non-trivial second-order approximate symmetry generators, 
but the first-order approximate symmetry generators are still retained. Thus there is no new 
approximate conservation laws but only the previous conservation laws that have been recovered. 

A better idea of what is actually required comes from the full system of geodesic equations. 
The geodesic equations are given by 

t + v'ir = 0, 

r + \{e v )'{e v c 2 i 2 - e~ u r 2 ) - re u {6 2 + sin 2 6>0 2 ) = 0, 

■■ 2 • -o 

6 + -r6- sinO cos 6»0 2 = 0, 
r 

0+ -r0 + 2cot(900 = 0, (16) 
r 
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with e 



now given by (|14p and hence 



v 



, e 1 - 2k 9 
= 1 e 

y 2 




c 



2k 



— i/ 



1 + - + 



e , (1-fe) 



2 



e 



r 



r 



r r 



where "•" denotes differentiation with respect to the geodetic parameter s. 

To construct the determining equations for the second-order approximate symmetries we 
use (|5|) and © as the 4 exact symmetry generators and ([T|)~ (fl~2]) as the 6 first-order approxi- 
mate symmetry generators. Of the 4 exact generators 2 do not appear in the new determining 
equations and the other 2 cancel out. The 6 generators of the first-order approximate symmetry 
have to be eliminated for consistency of the new determining equations, making them homoge- 
neous. The resulting system is the same as for Minkowski spacetime, yielding 10 second-order 
approximate symmetry generators. Four of them are again the exact symmetry generators used 
earlier, and hence simply "add into" them, making no difference. The other 6 replace the lost 
first-order approximate symmetry generators. The full set has the Poincare algebra so(l, 3)ffi s M 4 
apart from some non-Noether symmetries. There are no non-trivial second-order approximate 
symmetries as was the case for the first-order approximate symmetries. 

It is worth remarking that for the first-order approximate symmetries it did not matter 
whether we used the full system ([2]) or the un-perturbed system Eo, in ([4]). However, for the 
second-order approximate symmetries it does make a difference. One needs to use the full 
system ([2]) and not the un-perturbed system Eo, in (j4]) to obtain the solution. 

The exact symmetry generators include not only ([5]) and ©, but also the generators of the 
dilation algebra, d/ds, sd/ds corresponding to 

£{s) = c s + c\. 

In the determining equations for the first-order approximate symmetries the terms involving 
£ s = Co cancel out. However, for the second-order approximate symmetries the terms in £ s do not 
automatically cancel out but collect a scaling factor of (1 — 2k) so as to cancel out. (This factor 
comes from the application of the perturbed system, rather than the un-perturbed one.) Since 
energy conservation comes from time translational invariance and £ is the coefficient of d/ds in 
the point transformations given by ([I]), where s is the proper time, the scaling factor (1 — 2k) 
corresponds to a re-scaling of energy. Thus, whereas there was no energy re-scaling needed 
for the first-order approximate symmetries, it arises naturally in the second-order approximate 
symmetry. Using (|15[) we get the energy re-scaling factor 



Thus, even though there are no non-trivial second-order approximate symmetries, we get the 
non-trivial result of energy re-scaling from the second-order approximation. This point will be 
discussed further in the next section. 

It is worth remarking that when some symmetries are "lost" at one order (exact or first-order 
approximate) they are "recovered" at the next (at least to second-order) as "trivial" approximate 
symmetries. 

4 Summary and discussion 

We studied the approximate symmetries of the RN metric. This metric has isometry algebra 
so(3) ©M with generators ([5]) and ©. The symmetry algebra of the geodesic equations for this 
metric is so(3) © R © d^. We used Lie symmetry methods for DEs to explore the second-order 
approximate symmetries of the RN metric. Neglecting terms containing e 2 in the geodesic equa- 
tions (|16j) this metric has the same first approximate symmetries as those of the Schwarzschild 



(l-2fc) = (1 -Q 2 /2GM 2 ). 



Second-Order Approximate Symmetries of the Geodesic Equations 



7 



metric. Again we get no non-trivial approximate symmetry generator in the second approxi- 
mation. We only recover the "lost" conservation laws as approximate conservation laws. As for 
the Schwarzschild metric, where there is a difference between the conservation laws obtained for 
the system of geodesic equations and for the single orbital equation, the difference also holds 
for the RN metric. 

The re-scaling of energy for the RN metric, which does not appear for the Schwarzschild 
metric is of special interest. Notice that the pseudo-Newtonian formalism [391 140} l4"T | 142] gives 
re-scaling of force by (1 — Q 2 /rMc 2 ). The reduction is by the ratio of the electromagnetic 
potential energy at a distance r to the rest energy of the gravitational source. It is position 
dependent. That would not be reasonable for the energy in the field by itself. The scaling 
obtained here, (1 - Q 2 /2GM 2 ), is more reasonable as relating the electromagnetic self-energy 
to the gravitational self-energy, with the radial parameter, r, canceled out. 

It would be of interest to apply the definition of second-order approximate symmetries to 
the Kerr metric where there are only two isometries Xo = d/dt, X3 = d/d<j) and a non-trivial 
Killing tensor [2] to see if the same result of energy re-scaling holds. It would also be important 
to check if the energy in the gravitational field could be obtained by considering approximate 
symmetries of gravitational wave spacetimes. For this purpose one can first consider some non- 
flat static spacetime. Since gravitational wave spacetimes are non-static solutions of the vacuum 
Einstein field equations (EFEs), one should perturb the static spacetime with a time dependent 
small parameter and then look at the approximate symmetries of this gravitational wave-like 
spacetime. As gravitational wave spacetimes are solutions of vacuum EFEs, one would have 
to calculate the Weyl and stress-energy tensor for this gravitational wave-like spacetime (for 
which the stress-energy tensor is non-zero) to equal orders of approximation to see how much 
energy is contained in the gravitational field and how much in the matter field. To understand 
the actual problem one will then have consider an exact gravitational wave solution and look 
at its approximate symmetries. For the exact gravitational wave solutions a formula for the 
momentum imported to test particle is already available [43]. The comparison of the two results 
might enable one to identify a physically significant energy content in the gravitational wave 
spacetime. 
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